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1. Introduction 

There has been a great deal of interest and activity in percolation theory 
since the two Saint-Flour courses, [76, 111], of 1984 and 1996 reprinted in 
this new edition. We present here a summary of progress since the first 
publications of our lecture notes. 

The second edition of Percolation, [77], was published in 1999 as a fairly 
complete, contemporary account of 'classical' percolation (but with only lim- 
ited coverage of first-passage percolation). Books published since in the gen- 
eral area include [78, 79] and, for two-dimensional theory, [27]. There has 
been in addition a third Saint-Flour course on percolation, on the Wulff con- 
struction by Raphael Cerf in 2004, [40]. The 2004 review by Howard of 
first-passage percolation, [103], is useful. 

The principal targets of this short introduction are to indicate areas of sig- 
nificance that have received attention in recent years, and to include a useful 
bibliography. It is not intended as a comprehensive survey of the vast amount 
of material related to percolation and published since 1984/1996. We restrict 
our attention mostly to problems associated with finite-dimensional lattices, 
and we take quite a broad view of what constitutes a percolation or first- 
passage percolation problem. We offer apologies to the many authors whose 
work does not receive explicit mention here, and we trust sufficient pointers 
are present to guide the reader through recent developments in percolation 
theory. 

2. Aspects of first-passage percolation 

2.1. Limit theorems. The main 'ergodic theorems' for first-passage per- 
colation were established in their original forms by Hammersley and Welsh 
in 1965, and in more refined forms by others later; see [111] for the history 
prior to 1984. It is an important open problem to identify and prove the 
corresponding distributional limits. Only a limited amount is known beyond 
the Gaussian limit theorem of Kesten and Zhang for a critical model in two 
dimensions, [115]. They were able to exploit the special structure of dual 
cycles in critical percolation. 

One expects a Gaussian limit for passage-times along a long tube, but how 
wide can such a tube be? In recent work [46], Chatterjee and Dey have shown 
a Gaussian limit for d-dimensional tubes with width growing as a power of 
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length (see also [1]). There is evidence that the largest acceptable power is 
connected to the so-called fluctuation exponent (see Section 2.2). 

The weak limit is not expected to be Gaussian for unrestricted first-passage 
problems (see the corresponding discussion of Section 2.5 for last-passage 
percolation). Estimates of variance are first steps towards identification of 
distributional limits. Following earlier results by others, Benjamini, Kalai, 
and Schramm, [18], showed that the variance grows no faster than n/logn 
in a special first-passage process. Their methods and conclusions have been 
extended by Benaim and Rossignol, [16]. 

Large-deviations for passage-times have been discussed since [111] by Chow 
and Zhang, [50], and Cranston, Gauthier, and Mountford, [54]. Correspond- 
ing results for maximum flows are summarized in Section 2.4. 

LaGatta and Wehr have initiated a theory of 'Ricmannian' first-passage 
percolation. The shape theorem of [118] has been followed by a more detailed 
study of geodesies, [119, 120]. 

2.2. Fluctuation and wandering exponents. Some aspects of the geom- 
etry of first-passage percolation are not yet well understood. For example, 
how great are the fluctuations in fastest paths, when do asymptotic shapes 
have positive curvature, how many infinite geodesies exist? These intertwined 
questions have attracted a great deal of attention in recent years. We sketch 
certain recent advances concerning a fluctuation theory for fastest paths. 

For two vertices x, y, write T(x, y) for the passage-time from x to y, and 
.D(x, y) for the maximum deviation between the fastest path from x to y and 
the straight line-segment from x to y. For clarity, we assume that the passage- 
times have a continuous distribution, so that the fastest path between two 
given vertices is almost surely unique. 

Roughly speaking, the fluctuation exponent is the number x such that 
T(x, y) deviates from its mean by order \x — y\ x . The wandering exponent is 
the number £ such that D(x,y) is typically of order \x — y\^ . The universal 
relation x = 2£ — 1 is believed to hold, and furthermore x = §> £ = § in two 
dimensions. This relation is sometimes referred to as the KPZ relation, after 
the authors of [110] (see also [117]). 

Several possible definitions of the exponents have been discussed in the 
literature. Following earlier progress by others including Newman and Piza, 
[130], Chatterjee has succeeded in proving a version of the KPZ relation, 
under suitable definitions and subject to an assumption on the passage-time 
distribution, see [45]. His conclusions have been refined and improved by 
Auffingcr and Damron, [6]. 

2.3. Geodesies. A geodesic is a semi- or doubly-infinite path tt such that, 
for all vertices u,v € n, the section of n from u to v is shortest in the metric 
induced by the passage-times. Motivated by work with Stein on spin glasses, 
Newman asked in [128] for bounds on the number of geodesies that can exist 
in a d-dimensional first-passage process. He has conjectured that, subject 
to natural conditions, no doubly-infinite geodesies ('bi-geodesics') exist, and 
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there are infinitely many semi-infinite geodesies ('uni-geodesics'). Existence 
of bi-geodesics on Z 2 is equivalent to existence of non-constant ground states 
for the related disordered fcrromagnct, [129]. 

The existence (or not) of geodesies is bound up with the curvature of 
the boundary of the limiting shape of the growth process, and also with the 
question of coexistence of populations in a competition model. Results to 
date are incomplete. Interested readers are directed to the recent paper of 
Auffinger and Damron, [5], and to [25, 69, 98, 121, 124, 147, 148] and the 
references therein. 

The geometry of geodesies in planar Riemannian first-passage percolation 
is the subject of a recent project of LaGatta and Wehr, [119, 120]. 

2.4. The maximum flow problem. In two dimensions, the dual of first- 
passage percolation is the problem of finding the maximum flow in a certain 
(dual) random network. The maximum-flow problem is well posed in higher 
dimensions, and is dual to the problem of finding the minimum size of certain 
'cut-surfaces'. The last problem is much harder than the more usual first- 
passage problem, since tricky geometrical issues arise in stitching surfaces 
together. 

Partial progress was made by Kesten, [112], as reported in [111]. A number 
of authors have worked on the problem since, with substantial results by 
Zhang, [154, 155], Rossignol and Theret, [133, 134], and Theret, [145, 146]. 
Remarkable progress has been made over a period of years, culminating in a 
recent series of papers by Cerf and Theret, [42, 43, 44]. They consider the 
general situation of maximal flows in a scaling limit of n^ 1 !/ 1 as n — > oo. 
In the special case of a cube in d dimensions, these results yield that the 
maximum flow between opposite faces of the n-cube is asymptotic to <f>n d ~ l 
for some <\> = 4>(F)i where F is the distribution of a typical edge-capacity. 
In addition, the authors of these papers have proved lower and upper large- 
deviation theorems (of surface order n d ~ 1 , and volume order n d , respectively), 
and also that <f> > if and only if the atom F(0) is sufficiently small that the 
edges with weight do not percolate. 

2.5. Last-passage percolation. Johansson, [106, 109], showed the presence 
of the Tracy-Widom distribution within a certain percolation-type problem. 
While this is a last- rather than a /?rs£-passage problem, it is included here 
both for its intrinsic interest, and for its potential implications for weak 
limits of /irst-passage times. To each edge e of the square lattice Z 2 we 
allocate a random variable T e with the geometric distribution, P(T = k) = 
(1 — q)q , k > 0. Let G m ^ n be the suprcmum of the passage-times of paths 
from (1, 1) to (m,n) each of whose steps are either upwards or rightwards, 
and let 7 > 1. It is shown in [106] (see also [107]) that there exist functions 
/i(7, q), (7(7,(7), with known closed forms, such that (G 7 „. n — ^n) / (an 1 / 3 ) is 
distributed asymptotically as the Tracy-Widom distribution. The exponent 
^ occurs here just as in the problem of the longest increasing subsequence, and 
this is the conjectured value of the fluctuation exponent x m t w0 dimensions. 
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Johansson's work develops several beautiful connections: to so-called de- 
terminantal processes, [108], to the Kardar-Parisi-Zhang (KPZ) universality 
class, [52, 108, 110, 117], and to the 'totally asymmetric simple exclusion 
process' TASEP, [15]. 

Graham, [74], has explained how to use the theory of concentration to 
obtain the sublinearity of variance for last-passage percolation. We mention 
also the recent papers of: Hambly and Martin, [95], who consider edge- 
weights with heavy tails; Lin and Seppalainen, [125], concerning limiting 
shape; and Kcsten and Sidoravicius, [113], concerning directed last-passage 
percolation. 

2.6. Greedy lattice animals. To the vertices of a lattice we allocate ran- 
dom real-valued scores. A greedy lattice animal of size n is a connected 
subgraph of n vertices including the origin such that its aggregate score is 
maximal across this class. This variant of first-passage percolation was ini- 
tiated by Cox, Griffin, Gandolfi and Kesten, [53, 64]. A connection with 
Euclidean first-passage percolation was explored by Howard and Newman, 
[104]. The geometry of greedy animals has been investigated by Hammond, 
[96], where further references may be found. 

3. Percolation and disordered systems 

3.1. Conformality in two dimensions. Possibly the most prominent re- 
cent advance in probability theory is Schramm's construction of the family of 
random curves now termed Schramm— Lowner curves (SLE); see [135, 136]. 
The theory of SLE has revolutionized the mathematics of planar phase tran- 
sitions, and is a pillar of the bridge between probability and conformal field 
theory. 

In applying SLE to percolation, a spatial limit theorem is required in or- 
der to transform the critical process on a discrete grid into a process on 
a continuum. This limit theorem has been proved in essentially only one 
case, namely site percolation on the triangular lattice. A key ingredient for 
percolation was the so-called Cardy formula of [39] , which is based on the ar- 
guments of conformal field theory, and asserts a formula for the probabilities 
of open crossings of large domains in critical models. Cardy's formula was 
proved in 2001 by Smirnov, [140, 141], for the site model on the triangular 
lattice. The last paper included an outline proof of what has been called 
'full conformal invariance', and this has been fulfilled and amplified since by 
Camia and Newman, [35, 36]. 

It is a significant open problem to prove conformal invariance for an ex- 
tensive family of critical two-dimensional percolation processes. 

The theory of SLE has been described in the Saint-Flour notes of Wendclin 
Werner, [149], as well as in Greg Lawler's book [122], the recent surveys 
[11, 12, 123, 144], and the references therein. Proofs and elaborations of 
Cardy's formula have appeared in [27, 31, 79, 139, 150]. 
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Camia, Fontes, and Newman have asked about scaling limits for near- 
critical percolation, [32, 33]. A series of papers on this topic has been 
promised by Garban, Pete, and Schramm, of which the first is [67]. 

In contrast to the case of two dimensions, only little has been achieved since 
[76, 77] for critical percolation in (slightly) higher dimensions. For example, 
it remains a tantalizing problem to show the continuity of the percolation 
probability at the critical point. Similarly, there has been essentially no 
progress on the percolation phase transition in d dimensions with 3 < d < 18. 
When d > 19, we are in the domain of the lace expansion, see [77, Sect. 10.3]. 

3.2. Criticality and universality in two dimensions. A core technique 
for the study of near-critical percolation in two dimensions is the theory of 
box-crossings: when do they exist, and how may they be combined? The 
methods of Russo and Seymour- Welsh (RSW) have proved very useful in a 
variety of settings, including the proof of Cardy's formula and conformality. 
More recently, they have been used to prove aspects of universality for bond 
percolation on isoradial graphs. This large family of models includes but is 
much broader than the more familiar inhomogeneous square, triangular and 
hexagonal models. See [87, 88, 89] and the overview [80, Sect. 3]. 

RSW theory has become streamlined by the recent introduction of the the- 
ory of influence and sharp threshold, reviewed in [79, Sect. 4.5]. Examples, 
in addition to the universality results above, include the proofs by Bollobas 
and Riordan that the critical probability of site percolation on the Voronoi 
graph of a Poisson process is ^, [26, 28], and their exploration of the rela- 
tionship between self-duality and criticality for certain polygon models in two 
dimensions, [29]. 

Mention is made here of the recent proof by Schramm and Smirnov, [137], 
that any scaling limit of a planar percolation is a so-called 'black noise'. 

The connective constant of a lattice C is the number k such that the num- 
ber of n-step self-avoiding walks on C grows in the manner of K n . Identifying 
the value of k is akin to calculating the value of a percolation critical point. 
Nienhuis's long-standing conjecture that K = y/2 + \/2 for the hexagonal lat- 
tice has been proved in notable work of Duminil-Copin and Smirnov, [61]. 
It is an important open problem to show that the weak limit of a random 
self- avoiding walk in two dimensions is the Schramm-Lowner curve SLE 8 / 3 
(see [10]). 

3.3. Dynamical percolation and noise sensitivity. Suppose a critical 
percolation model is allowed to evolve in time: the state of each bond/site 
switches between open and closed after exponential holding times, with rates 
chosen in such a way that the process remains critical. This 'dynamical per- 
colation' process was introduced in independent work of Haggstrom, Peres, 
and Stcif, [94], and of Benjamini. The main issue has been to determine the 
Hausdorff dimensions of the sets of times at which certain forms of aberrant 
behaviour occur. It has been normal to consider the case of two dimensions, 
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and particularly (but not exclusively) site percolation on the triangular lat- 
tice, since a certain amount is then known about pivotal sites for long open 
connections. Earlier results of Schramm and Steif, [138], for the existence of 
an infinite open cluster have been complemented by a more detailed analy- 
sis by Garban, Pete, and Schramm, [66]. Dynamical percolation has been 
reviewed recently by Garban and Steif, [68]; see also [30, 65, 143]. 

In the related area of 'noise-sensitivity' of Boolean functions, Benjamini, 
Kalai, and Schramm have asked how sensitive is the probability of an event 
to the resampling of the states of bonds/sites, [17]. Noise sensitivity has been 
reviewed by Garban, [65], and jointly with Steif, [68]. 

3.4. Incipient infinite cluster. As remarked in Section 3.1, it remains an 
open problem to prove that there is no infinite cluster in the critical perco- 
lation model in a general number of dimensions. This is known when d = 2 
(by box-crossing arguments) and when d > 19 (using the lace expansion). 
The so-called 'incipient infinite cluster' (IIC) is a measure defined via con- 
ditioning at or near the critical point, under which an infinite cluster exists. 
Kesten's proof of existence of the IIC in two dimensions has been extended 
and developed by Jarai, [105]. Corresponding results in high dimensions have 
been established by van der Hofstad and Jarai, [100], using the lace expan- 
sion, and more recently by Heydenreich, van der Hofstad, and Hulshof, [97], 
assuming a form of the triangle condition. The latter results are valid also 
for sufficiently spread-out models in 7 and more dimensions. 

Much is known about the IIC for oriented percolation, see Section 3.8. 

3.5. Random walk on the infinite cluster. Supercritical percolation pos- 
sesses an almost surely unique infinite cluster. This cluster has many proper- 
ties in common with the entire lattice, and these two sets may be said to be 
'coarsely' equivalent. One test for coarse equivalence is whether or not ran- 
dom walk on the infinite cluster satisfies an invariance principle. A number 
of authors have worked on this problem, with notable success. No complete 
bibliography is attempted here. We refer the reader instead to the survey by 
Biskup, [24], of the more general topic of random walks on lattices endowed 
with random conductances. 

There is an old conjecture of Alexander and Orbach for the rate of decay 
of the return probability of random walk on the incipient infinite cluster, 
namely that P2n(x, x) is of order n -2 / 3 . This has been proved by Kozma and 
Nachmias, [116], for large d (d > 19 suffices) and for the spread-out model 
with d > 6. A related result for oriented percolation is found in Section 3.8. 

We mention also recent work of Fribergh and Hammond, [62]. In a contin- 
uation of earlier work by others, they have studied the velocity of a random 
walk with drift on the supercritical infinite open cluster. They have proved 
the (arguably counter-intuitive) existence of a critical value for the drift rate 
below which the velocity is strictly positive, and above which it is 0. 



PERCOLATION SINCE SAINT-FLOUR 



7 



3.6. Subcritical percolation. In the subcritical phase, the two-point con- 
nectivity function decays exponentially with distance n. That is, the proba- 
bility of an open connection between two points at distance n behaves quali- 
tatively in the manner of e~ n ^^ where is the 'correlation length'. There 
is a multiplicative correction term n^^ 1 ^ 2 , and the ensuing decay is termed 
'Ornstein-Zernike decay' (OZ), [131]. This polynomial correction has been 
proved throughout the subcritical phase by Campanino and Ioffe, [38], in a 
refinement of earlier work restricted to axial directions. In the supercritical 
phase, one considers instead the truncated two-point function, for which the 
most recent study is that of Campanino and Gianfelice, [37]. It is an open 
question to prove OZ decay for the truncated function all the way down to 
the critical point. 

3.7. Supercritical percolation. A finite open cluster in two dimensions is 
surrounded by a dual cycle. Duality is much more complex in dimensions 
exceeding 2, since the boundary of a finite cluster is a surface, and general 
surfaces have complicated geometry The geometry of random surfaces has 
been the source of a number of fine problems, of which we mention three. 

Whereas classical percolation concerns the binary relation of 'connectiv- 
ity', one may instead seek conditions for there to be an infinite 'entangled' 
set of open edges (in three dimensions). Unlike 'connectivity', entanglement 
is a long-range property, and its study involves novel complications of a ge- 
ometrical nature. Early papers [82, 92, 101] by Grimmett, Haggstrom, and 
Holroyd stimulated recent work of Atapour and Madras, [4], who have proved 
an exponential-tail theorem for the size of the largest entanglement at the 
origin in part of the subcritical phase. An improved lower bound for the 
critical entanglement probability has been proved by Grimmett and Holroyd, 
[83]. 

The last reference is the second of a sequence of four papers by the same 
authors, sometimes joined by Dirr, Dondl, and Schcutzow, devoted to ques- 
tions associated with so-called 'Lipschitz percolation', see [57, 84, 85]. Lips- 
chitz percolation is concerned with the existence and geometry of 'flat' open 
surfaces within a percolation process. 

One may ask similarly for conditions under which there exists an infi- 
nite 'rigid' set of open edges. This problem has been studied by Holroyd, 
[102], and Haggstrom, [91]. In [51, 127], Connelly, Menshikov, Rybnikov, 
and Volkov have proved loss of tension in planar networks subjected to ran- 
dom perforations. 

By conditioning on the absence of open connections between the upper 
and lower hemispheres of a large sphere, Giclis and Grimmett, [72], have con- 
structed Dobrushin interfaces for percolation in three and more dimensions, 
and more generally for random-cluster models. The associated 'roughening 
transition' is no better understood for such systems than for the Ising model. 

The chemical distance between two points is defined to be the length of 
the shortest open path between these points. Antal and Pisztora showed 
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the 'law of large numbers' for chemical distance, [3]. Their work has been 
elaborated recently by Garet and Marchand in their studies of large and 
moderate deviations, [70, 71]. 

3.8. Oriented percolation. The existence of the incipient infinite cluster 
for oriented percolation has been proved in d > 4 + 1 dimensions by van der 
Hofstad, den Hollander, and Sladc, [99]. 

In a useful paper concerning random walk on random graphs, Barlow, 
Jarai, Kumagai, and Sladc, [8], have proved a version of the Alexander- 
Orbach conjecture in the setting of the incipient infinite cluster of oriented 
percolation in d > 7 dimensions. 

The geometrical structure of supercritical oriented percolation has been 
investigated by Grimmett and Hiemer, [81], and further in two dimensions 
by Wu and Zhang, [152]. 

3.9. Invasion and bootstrap percolation. 'Invasion percolation' is a sys- 
tem that is "self-organized" to focus on the critical point of percolation. It 
proceeds by a progressive filling of so-called ponds. Quite a lot is now known 
about these ponds and their outlets. The reader is referred to the recent pa- 
pers [55, 56] by Damron and Sapozhnikov, whose results have been inspired 
in part by earlier work of others, [2, 73], on invasion percolation on trees. 

In 'bootstrap percolation', each vertex in a box of side- length n is declared 
'infected' with probability p, and further vertices are infected if they have two 
or more infected neighbours. The critical point p c = p c (n) is the value of p 
above which there is probability at least \ that every vertex becomes infected. 
Increasingly accurate asymptotics of p c {n), viewed as a function of n as n — > 
co, have been obtained. Gravner, Holroyd, and Morris, [75], have provided 
the second-order term in the expansion of p c {n) in two dimensions, thereby 
complementing the earlier result of Holroyd that p c (n) ~ 7r 2 /(18 logn). The 
sharpness of the phase transition for a general class of such models has been 
established by Balogh, Bollobas, Duminil-Copin, and Morris, [7]. 

3.10. Percolation of words. Each vertex of a lattice is allocated a random 
letter from the two-letter alphabet {0, 1}. A word is a semi-infinite sequence 
of letters. Benjamini and Kesten, [19], initiated the study of the set of words 
seen in a two-dimensional lattice. Their work has been continued by Kesten, 
Sidoravicius, and Zhang, [114]. 

Grimmett and Holroyd, [85], have presented related results for Lipschitz 
embeddings of <i-dimensional words in <i'-dimensional space. The case d = 
d! = 1 was posed by Grimmett, Liggett, and Richthammer, [86], and has 
been largely solved in recent work of Basu and Sly, and of Gacs, [9, 63]. 

3.11. Percolation on non-amenable graphs. In the last 20 years, there 
has been a good deal of interest in random processes on graphs whose sur- 
face/volume ratio does not tend to zero. Grimmett and Newman, [90], proved 
the existence of three phases for percolation on the direct product of a tree 
and a lattice. This provoked Benjamini and Schramm, [21], to propose in 
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1996 a systematic study of percolation models on non-amenable graphs. They 
made a number of conjectures, some of which have since been resolved. For 
a recent review of this field, the reader is referred to the memoir, [93], by 
Haggstrom on the contributions of Odcd Schramm. 

The nature of percolation on a graph G is connected to the algebraic 
structure of G and its automorphism group. Of particular interest are Cay- 
ley graphs of finitely-generated groups. The mass-transport principle was 
introduced to percolation for graphs with unimodular automorphism group 
by Benjamini, Lyons, Peres, and Schramm, [20]; see [93, Sect. 3]. 

3.12. Disjoint occurrence. The BK inequality is one of the main tools of 
percolation theory. In its simplest form, it assumes product measure and it 
applies to increasing events. The 'increasing' assumption was removed by 
Rcimcr, [132]. The assumption of independence has proved resilient, but 
has been weakened by van den Berg and Jonasson, [23], to randomly drawn 
subsets of given size. This in turn has provoked a further extension to certain 
Ising-type models by van den Berg and Gandolfi, [22]. 

3.13. Random-cluster and Potts/Ising models. Geometrical techniques 
have long been used in the study of phase transitions. The random-cluster 
representation is a prominent example, in which the two-point correlation 
functions of the Ising/Potts models are transformed into the connection prob- 
abilities of a dependent model of percolative type. Recent books on the 
random-cluster model include [78, 151]. 

There has been major progress on the theory of random-cluster models, 
although many questions remain open. One significant advance concerns the 
rigorous theory of the so-called Wulff construction. Wulff, [153], presented 
a generic formula for the shape of a large 'droplet' in a statistical physi- 
cal model, as the 'crystal' of given volume with smallest aggregate surface 
tension. Cerf and Pisztora, [41], have formulated and proved the Wulff con- 
struction for random-cluster models in general dimension. Their solution, 
together with Cerf 's solution of the Wulff crystal for percolation, is the sub- 
ject of Raphael Ccrf's Saint-Flour notes, [40]. 

Of the many open problems in this area, we mention one. No proof is yet 
known of the uniqueness of infinite- volume measures above the critical point 
in three or more dimensions, [78, Thm 5.33]. 

Just as the two-dimensional Ising model has proved very special, so are 
planar random-cluster models with cluster-weighting factor q = 2. Smirnov, 
[142], has shown conformal invariance of the critical model on the square 
lattice, and this work has been extended (jointly with Chelkak) to isoradial 
graphs, see [48, 49] and the lecture notes [60]. Duminil-Copin, Honglcr, 
and Nolin have extended RSW-type results to this model, [59], and this 
has provided a key technique for the proof by Lubetzky and Sly, [126], of a 
polynomial bound for the mixing time of the critical Ising model. It is still 
unclear to what degree techniques for q = 2 may be exploited for general 
q > 1. It has been shown in [14] that the so-called 'parafermionic observable' 
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of [142] may be used to study the q > 4 random-cluster model on the square 
lattice. 

In a further study of the critical Ising model on the re-scaled square lattice 
aZ 2 , Camia, Garban, and Newman, [34] , have shown that the magnetization 
field, scaled by a 15 / 8 , converges as a — > to a conformally covariant field 
which is non-Gaussian. Chelkak, Hongler, and Izuryov, [47], have proved the 
asymptotic conformal covariance of the n-point correlation functions of this 
critical Ising model. The near-critical FK-Ising and random-cluster models 
have been considered by Duminil-Copin, Garban, and Pete, [58]. 

The exact critical value of the random-cluster model on the square lattice 
has long been "known but unproven" . This has now been rectified by Beffara 
and Duminil-Copin, [13], in their proof that p c (q) = <Jq/(l + <Jq) for q > 1. 
Box-crossings play a central role in their arguments. See also [14, 80]. 
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